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We consider the Bose-Hubbard model with particle losses at one lattice site. For the non- 
interacting case, we find that half of the bosons of an initially homogeneous particle distribution, 
are not affected by dissipation that only acts on one lattice site in the center of the lattice. A 
physical interpretation of this result is that the surviving particles interfere destructively when they 
tunnel to the location of the dissipative defect and therefore never reach it. Furthermore we find 
for a one-dimensional model that a fraction of the particles can propagate across the dissipative 
defect even if the rate of tunneling between adjacent lattice sites is much slower than the loss rate 
at the defect. We analyze the robustness of our findings with respect to small interactions and small 
deviations from the symmetric setting. A possible experimental realization of our setup is provided 
by ultracold bosonic atoms in an optical lattice, where an electron beam on a single lattice site 
ionizes atoms that are then extracted by an electrostatic field. 
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I. INTRODUCTION 

Ultracold atoms in optical lattices have in recent years 
emerged as a very successful quantum simulator for 
many-particle and solid state physics. The ability to 
tune parameters of the simulated Hamiltonians provides 
a unique tool to explore quantum phase transitions in this 
system 1 . Among the most studied models is the Bose- 
Hubbard Hamiltonian which describes interacting bosons 
in a lattice potential that can tunnel between neighbor- 
ing lattice sites. In its ground state, it can show two 
quantum phases, a superfluid regime for weak interac- 
tions and a Mott-insulator phase for strong interactions 
and commensurate filling [2]. 

Electron beams can provide means for locally prob- 
ing ultracold quantum gasses since they can be focused 
onto spot sizes that are much smaller than the optical 
wavelengths of the trapping fields. High-resolution scan- 
ning electron microscopy thus constitutes a method that 
combines single atom sensitivity and high spatial resolu- 
tion H] . For this technique, a focused electron beam 
scans an optically trapped atomic gas. The atoms are 
ionized by the electron impact and the produced ions 
are moved off the dipole trap by an electrostatic field. 
With this method, single-site addressability of ultracold 
atomic gases loaded in optical lattices has been demon- 
strated in one and two dimensional structures [H |S] . In 
addition, second and higher order correlation functions 
of a trapped gas of bosons have been extracted [6] . 

Besides providing a tool for local measurements, a fo- 
cused electron beam can also be used as a method to 
introduce controlled localized dissipation. Along these 
lines, the formation of solitons in a continuous system 
with a localized dissipative perturbation has been dis- 
cussed employing a mean field approach [7] and correc- 
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tions beyond the mean field level have also been con- 
sidered [5]- Moreover, a three-site Bose-Hubbard model 
with dissipation in the central lattice site has been ana- 
lyzed [3J and a numerical study of the dynamics during 
a finite time range that employed DMRG methods has 
been presented [TP] . 

In this article we investigate the effect of a localized 
dissipative defect in a Bose-Hubbard model with analyt- 
ical means. This allows us to consider arbitrarily large 
lattices as well as the long time limit of the dynamics. For 
the non-interacting case, we show that a certain fraction 
of the population of a uniform particle distribution is not 
affected by dissipation that is localized in the center of 
the trap. In other words, we show the existence of a 
dissipation-free subspace for localized particle losses. 

Furthermore we find for a one-dimensional model that 
a fraction of the particles can propagate across the dissi- 
pative defect even if the rate of tunneling between adja- 
cent lattice sites is much slower than the loss rate at the 
defect. Particles can however only tunnel through the 
dissipative defect at the expense of other particles being 
lost via the localized dissipation. 

For the description of the system we employ a Born- 
Markov master equation in the standard Lindblad form 
[TTHH1- This equation represents the dynamics of the 
system under the influence of a mechanism that intro- 
duces one-by-one particle losses at one lattice site. 

We start by investigating a non-interacting system 
with the defect located exactly at the center of the chain. 
In this case, we show the existence of a dissipation-free 
subspace which consists of the half of the normal modes 
of the Hamiltonian. Those are the modes that correspond 
to wave-functions which are anti-symmetric with respect 
to the trap center. The bosons of those modes interfere 
destructively at the location of the defect and therefore 
never reach the lossy site. Since their wavefunction is fi- 
nite everywhere else, these bosons can tunnel between the 
right and the left sides of the system without appearing 
at the defect. 
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After this observation, we investigate cases that might 
present dynamical behavior similar to that of the non- 
interacting and symmetric one. Most important candi- 
dates are cases that correspond to small deviations from 
the non-interacting and symmetric case. One way to 
introduce such a deviation is to turn on small interac- 
tions while keeping the symmetric structure of the sys- 
tem. Interactions destroy the phase coherence along the 
lattice since they can transform protected to unprotected 
modes. Also in a non-interacting system with a slightly 
non-symmetric structure, we find that the destructive in- 
terference of particles that tunnel into the defect and 
therefore the dissipation free subspace are perturbed. 
For those two types of small deviations from the non- 
interacting and symmetric case, namely small interac- 
tions and slightly non-symmetric structure, we show that 
there is a slow down of the losses for some fraction of the 
total population. 

The remainder of the paper is organized as follows. 
In section II we introduce our model and discuss gen- 
eral properties of particles in trapping potentials that 
are symmetric around a dissipative defect. In section III, 
we present quantitative results for non-interacting bosons 
on a lattice with a dissipative defect in the central lat- 
tice site. Section IV then discusses small deviations from 
the interaction-free and perfectly symmetric case and we 
finish the presentation with a summary and conclusions 
in section V. 



II. LOCALIZED LOSSES IN SYMMETRICAL 
POTENTIALS 

Various quantum systems can feature states that are 
protected against certain types of dissipation. Those 
states - often referred to as dark states - are well known 
in atomic physics and quantum optics [TTJ fTS] . In many- 
body systems such states can appear as a property of a 
symmetry of a given Hamiltonian and be part of its eigen- 
states. One of the symmetries that can generate such a 
dissipation-free subspace is the symmetry under reflec- 
tion at the origin of the system. In other words, when 
the potential is symmetric around a localized dissipative 
defect. 



Dissipation-free subspace of a single particle in a 
symmetrical potential 

We consider a single particle in a symmetric poten- 
tial. Such a potential is described by an even function 
U with U(x) — U(—x). The wavefunctions that cor- 
respond to eigenstates of such a potential must have a 
well defined parity, i.e. some are even (ip e (x) — ipe{~x)) 
and others are odd (jp {x) — —ipo{—x)). These proper- 
ties emerge since the Schrodinger equation is unchanged 
when the sign of a coordinate is reversed |15j . If one lists 
the eigenfunctions according to the number of nodes - in- 



creasing energy - then they are alternately even and odd. 
The ground state of the system must be an even function 
since it contains no nodes. Odd functions always vanish 
at the origin, since 

i>o(Q) = -1>o(o) = o. (i) 

This means that whatever happens at this particular lo- 
cation in space does not affect the odd eigenfunctions. A 
physical interpretation of this result is that the probabil- 
ity amplitudes for particles in odd eigenstates destruc- 
tively interfere at the origin. 

Now, let us assume that we locate a mechanism that 
introduces particles losses at the origin. Then, if the par- 
ticle is in a odd eigenf unction, it will never be caught by 
the loss mechanism. This is to say, the odd eigenfunc- 
tions constitute a dissipation-free subspace. 

Extension to a lattice setup 

The effect of the dissipation-free subspace may seem 
(almost) unphysical in a continuous description. This 
is because we silently assumed that the loss mechanism 
has zero spatial width. However, this wouldn't be prob- 
lematic in a lattice setup. If we assume that the rel- 
evant physics takes place within the lowest band of a 
tight binding Hamiltonian, i.e. the energy of the sys- 
tem is low enough, then the lattice sites can be treated 
as points in space since their extension is much smaller 
than the wavelength of the mode-functions. Therefore, 
the only requirement for the loss mechanism is to be re- 
stricted to one lattice site, specifically the central one 
- we assume odd number of sites. Therefore the mech- 
anism that generates the losses needs to be capable of 
addressing single lattice sites, see Fig. (JTJ. Such mecha- 
nisms have recently been realized for ultracold atoms in 
optical lattices [JJ [H HB], which motivates the study 
we present here. 

Description of the losses 

We assume that particle losses occur one-by-one and 
can be described by a super-operator within the Born- 
Markov approximation [Til H2] which reads, 

\ps(t)] = x (2a>nPs{t)al - ala n ps(t) - ps(t)a^a n ) , 

where ps(t) is the (reduced) density matrix of the sys- 
tem and ajj and a n are the creation and annihilation op- 
erators of one atom at site n. For experimental realiza- 
tions of local particle losses induced by a focused electron 
beam, our assumptions for the decay term are well justi- 
fied. Here the atoms are ionized by the electron impact 
and then extracted from the lattice with an electric field. 
Therefore subsequent particle decays will be independent 
of each other which justifies the Markov assumption, c.f. 
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FIG. 1: (Color online) Illustration of the proposed setup for 
the realization of the dissipation-free subspace. Introducing 
localized particle losses at the central site - assuming odd total 
number of lattice sites - leaves the odd modes - described by 
anti-symmetric functions - unaffected. The loss mechanism 
consists of an electron beam focused onto the central site. 
The atom-electron collisions leave the atoms ionized and by 
the use of an electrostatic field they can be extracted from 
the lattice. This method provides single-site sensitivity for 
the control over the position of the defect and high precision 
control of the associated damping rate. 



wise) on-site particle interaction. The operators &i and 
a\ are the bosonic annihilation and creation operators, 
respectively. They satisfy the commutation relations for 
bosons, [aj,oJ] = Sij and [04,07] = [oJ,o]-] = 0. Since we 
are interested in the dynamics generated by the Hamil- 
tonian (|3| for a fixed number of particles, we omit a term 
related to a chemical potential. Two quantum phases of 
this model are well known, namely the superfluid and the 
Mott insulator [5] . In general the phase of the sample de- 
pends on the magnitude of the ratio J/U. An experimen- 
tal realization of such a model is provided by ultracold 
bosonic atoms in an optical lattice [TJ . 

Feshbach resonances provide a unique and efficient way 
to control interactions. Based on this method, an ideal 
system of non-interacting atoms can be approximately 
implemented. In this case Eq. ^ is approximately writ- 
ten 

H = -jJ2(4ai + h.c.y (4) 

(id) 

This is the hopping Hamiltonian, and because of the ab- 
sence of interactions the dynamics of each particle is un- 
affected by presence of the remaining particles. 



Addition of a confining potential 



A. The dynamics in one dimension 



In an experimental setup with cold atoms in an opti- 
cal lattice, the atoms are not only subject to the lattice 
potential but are also trapped by an additional confining 
potential that is usually well approximated by a harmonic 
trap. In the central region of the lattice, the value of the 
confining potential typically only varies very little from 
one lattice site to another and we therefore approximate 
it by a constant. In current experiments it is however 
not guaranteed that the bottom of the confining poten- 
tial exactly overlaps with the site of the lattice where 
losses are induced. To estimate the effect of this feature 
on our findings, we consider deviations from a lattice that 
is symmetric with respect to the dissipative site in section 

WJ 



III. NON-INTERACTING BOSONS IN A 
LATTICE WITH PARTICLE LOSSES AT THE 
CENTRAL SITE 



For simplicity, we start with the dynamics in one di- 
mension. However, as we shall see, the generalization 
to higher dimensions is straight forward. The Hamilto- 
nian Q in one dimension is written in the form 



JV-l 



H = - J (a^a„ + i + h. 



(5) 



where N is the total number of lattice sites. We assume 
open boundary conditions, which imply that the wave- 
function vanishes at the edges of the system. The ladder 
operators for the normal modes, bk, are in this case given 
by the following transform 



N- 



(6) 



N ■ 



— ^ sin (nk)a n , 



We start our quantitative study by considering a non- 
interacting gas of bosonic particles on a lattice with par- 
ticle losses at the central site. Assuming that the tight 
binding approximation holds, the model under consider- 
ation is the Bose-Hubbard Hamiltonian, 



H = -J 



(id) 



U 
2" 



^ota* - i) > ( 3 ) 



where J is the tunneling rate of the bosons between 
neighboring sites and U is the strength of the (pair- 



where the quasi-momcntum k is defined as k — , with 

I = 1, 2, , N. It can be easily shown that the operators 

bk and frj. satisfy the commutation relations for bosons. 
By substituting Eq. Q into Eq. Q, we can write the 
Hamiltonian in the diagonal form 



H 



-2J^2cos(k)blb k . 

k 



(7) 



This Hamiltonian corresponds to a set of N harmonic os- 
cillators with frequencies given by the dispersion relation 
uj(k) = — 2Jcos(k). 
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We note that even though the momentum modes ^ 
only diagonalize the Hamiltonian ^ in the absence of 
a confining potential, there is a different set of normal 
modes which diagonalizes the model with harmonic con- 
fining potential as there arc no interactions. Following 
our discussion in section |n] half of those normal modes 
will nonetheless be protected against the localized dissi- 
pation. 



The Master Equation 

The Born-Markov master equation that describes the 
evolution of the system under the influence of the dissi- 
pation at one lattice site is of the form, 



dt 



Ps(t) 



H,p s (t) +C m [ Ps (t)} 



(8) 



where we have set H = 1. The dissipator in the above 
master equation describes particle losses at site to, as 
defined in Eq ([2]). In a first approximation, we shall 
assume that this site corresponds to the central site of 
the lattice. This assumption of a symmetric structure is 
motivated by the fact that the lattice site on which the 
electron beam is focused can in experiments be selected 
with great accuracy [H [5] . 

In terms of the mode operators given by Eq. the 
dissipator of the master equation reads 



£m [f>s{t)} 



N ■ 



— ^ sin ( mfc ) sin ( m P) { 2 i>kPs(t)bl 



(9) 



Under the assumption that the total number of sites is 
odd, the central site is numbered as to = [N + l)/2 
and by using the definition of the quasi-momentum, the 
dissipator in Eq. ^ is written 



£m[p S (t)}=^ L 1 j: 

Lj = i 

b}bjPs(t) - Ps(t)b]bj} 



N 



sm l y) sin (¥) {2/ ' ,/M/i/ ' 



(10) 



We immediately see that only odd values of I and j con- 
tribute to the summation. That means that the dissipa- 
tion affects only the modes that correspond to odd quasi- 
momentum labeling number. Therefore, the modes with 
even constitute a dissipation- free subspace of the dis- 
sipator as given in Eq. 

At first site, this result seems to be in contradiction 
with the discussion in the previous section. There, we 
concluded that the dissipation-free subspace consists of 
the odd eigenfunctions. However, this contradiction is a 
misconception since the odd eigenfunctions correspond to 
even quasi-momentum labeling number and vice versa. 



The master equation, given by Eq. ^ can now be 
written in the following form, 



dt 



Ps(t) 



N 



H,p s {t) 

N 



7 



mii I ? Wn(^ ) {2bip s {t)b 



b\hps{t)-Ps{t)b]b 3 }. 



(11) 



Using the above equation, one can derive equations 
governing the evolution of the mean values of observ- 
ables. For the mean value of the number of bosons in 
odd modes (with even quasi-momentum labeling num- 

ber) (N(t)) 0d d = EAblMit) = Ej' Tr{bl v b 2l ,p s (t)}, 
where I' = 1, 2, N/2, one gets, 



dt 



(N(t)) 



Odd 



0. 



(12) 



The obvious solution of this equation, (N(t))odd — 
(NlO^Qdd confirms that modes with an even quasi- 
momentum labeling number constitute a dissipation-free 
subspace. 

On the other hand, for the mean value of the num- 
ber of bosons in even modes (with odd quasi-momentum 
labeling number) (N(t)) Even = J2i> (^L'-i^'-i)^) = 

Ei' Tribly^kv-iPsit)}, where I' = 1, 2, (N + l)/2, 
we have, 



dt 



(N(t)) 



Even — 



27 
N+ 1 



(N+l)/2 



(*) 



(13) 



Z,j=l 



Hence, particles in these modes will eventually be lost 
from the system. The numerical solution of the above 
equation is shown in Fig. One sees that all bosons 
that populate the even modes eventually decay. 



Particle propagation across the defect 

We now consider the case where initially one boson is 
located at a specific lattice site at one side of the dissi- 
pative defect. The state of this localized particle is a co- 
herent superposition of all the eigenmodes of the system 
[17] , Therefore, the fraction of the particle that occupies 
odd modes will be gradually lost due to the dissipation at 
the defect whereas the remaining part will be protected 
against the localized losses as it occupies even modes. 
Since the even modes however show equal particle densi- 
ties at both sides of the defect, a fraction of the initially 
localized particle is able to propagate through the defect 
into the other half of the lattice, see Fig. In complete 
analogy, a fraction of a gas of non-interacting bosons that 
is initially located at one side of the dissipative defect, 
will propagate through the lossy site into the other half of 
the lattice at the expense of loosing half its initial number 
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FIG. 2: (Color online) The expectation value of the popula- 
tion of the even modes as a function of time. All the bosons 
that populate these modes eventually decay. The plot reveals 
also the rapid loss of the particles initially populating the de- 
fect. A large proportion of the particles is annihilated almost 
instantly in the case where 7 ^> J. This suggests that, in this 
case, the lossy site can be adiabatically eliminated from the 
description because of its fast rotating dynamics. Parameters: 
7 = 10, J = 0.1, U = and N = 5. 



B. Extension to higher dimensions 

Since in our case the dissipation- free subspace is a con- 
sequence of the symmetry of our system, one expects to 
find an analogous behavior in higher dimensions. To this 
end, we consider a two dimensional square lattice of size 
N x N and suppose that the electron beam focuses ex- 
actly on the center of the square lattice. In this case, the 
Hamiltonian given by Eq. Q takes on the form, 



JV-l 



H 



E (4: 



+ at Aij+i + h.c. 



(14) 



By assuming again open boundary conditions, the ladder 
operators are transformed as, 



1 -3 



N ■ 



— ^2 sin ( ik ) sin Up) i>k, P , 



k , p 
N 



(15) 



bk, P — 



N ■ 



— ^2 sin (ik) sin (jp)a itj , 
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where again the quasi-momenta take the values k,p = 
j^ri, with I = 1, 2, N. In terms of the mode operators 



bk, P and b\ p , Hamiltonian (j 1 4|) has the following diagonal 
form, 



H = -2JJ2 I cos ( k ) + cos (P)l b iJk, P - 

k . p 



(16) 



The dissipator that describes the particle losses at the 
central site reads, 



7 . 

£m,m [PS\t)] =2 [2<im,mPs(t)a. 



t 

m.m 



Ps(t) 



m,m a m,mj 



(17) 



FIG. 3: (Color online) The propagation across the defect 
of a particle that is initially localized at the first lattice site 
from the left. Dashed line: The expectation value of the total 
population in lattice sites to the left of the defect. Thin line: 
The expectation value of the total population in lattice sites 
to the right of the defect. Thick line: The expectation value 
of the population in the defect. Half of the initial number of 
particles survives the damping and ends up in a superposition 
of being located on the left and the right part of the lattice. 
Parameters: 7 = 1, J = 0.1, U — and N = 5. 



of particles, see Fig. Q. It is remarkable that this effect 
is independent of the rate of particle losses at the defect. 
Hence even if the particles are lost from the defect a lot 
faster than they can tunnel into it, still a fraction of the 
particles (1/4 of them in our examples) is able to cross 
the defect and reach the other half of the lattice. 



where the coordinates of the lossy site are given by m — 
(N + l )/2. In terms of the eigenmode operators defined 
in Eq. (151, the above dissipator reads, 



IPs®] 



27 



(N+ l) 2 

N 



• I m \ ■ ( n 3 \ ■ ( nl \ ■ I n f 
> sin I — I sin I — ] sin | — | sm 



x{2b. hJ p s (t)bj f - b}jbi,fP8(t) - PsitWiJlj}- 

(18) 

One immediately sees that the above dissipator is non- 
zero only when both the entries of the summations 
are odd (corresponding to even modes). The eigen- 
modes that correspond to any other combination of 
quasi-momentum labeling numbers, constitute part of 
the dissipation-free subspace. This can be demonstrated 



G 



by computing the evolution of the total number of par- 
ticles in modes corresponding to every such combina- 
tion. For the modes where at least one index of quasi- 
momctum labeling is an even number, one obtains, 



j t {N(t))o-o = ~{N(t)) E -o 



0. 



dt 



(N(t)) 



E-E 



4-y 



{N + iy 



(N+l)/2 
l,j,kj=l 



tween the defect and its neighbor sites, 



V(t) = — J [a m _ 1 (t)a m (t) + a] n+1 {t)a m {t) + h.c. 



We now focus on the regime 7 > J where we can as- 
sume that the lossy site practically remains empty at all 
times, as illustrated in Fig. ( 2|). By imposing the Born 
approximation the equation (211) reads, 



(19) 



On the other hand for the modes where both indices are 
odd numbers (even modes) one gets, 



l,2i-l°2fc-l,2/- 17(f) 

(20) 

We therefore conclude that the dissipation-free subspace 
of a two dimensional non-interacting, symmetric system, 
constitutes the 3/4 of all its cigenmodes. 

In a completely analogous way, one finds that in a 
three-dimensional structure, the protected modes are 7/8 
of all its eigenmodes. 



The case of large damping rate and the 
adiabatic elimination of the defect 



Let us return to the one dimensional setup. As illus- 
trated in Fig. ([2]), in the case where the damping rate 7 
is much larger than the tunneling J, the loss of the par- 
ticles that populate the defect, occurs almost instantly 
compared to the time-scale in which the tunneling takes 
place. This suggests that one can assume the lossy site 
to be empty at all times and adiabatically eliminate its 
degrees of freedom to generate an effective picture for 
the description of the remainder of the system. We now 
discuss this approach. 



The effective Markovian dissipator 

We consider a system with three sites [9] with particle 
losses occurring at the central site. By performing per- 
turbative theory in the interaction picture |111 112j and 
tracing out the degrees of freedom of the lossy site, we 
derive an expression for the evolution of the reduced den- 
sity matrix of two outer sites p_ m , 



C eff [p- m (t)} = - Tr m V(t), V(t-s),p(s) ds. 

Jo 1 1 JJ 

(21) 

Since we trace out the degrees of freedom of the defect, 
we automatically place the lossy site in the status of a 
quantum environment. Thus, the interaction Hamilto- 
nian V(t) in Eq. (21 1 corresponds to the tunneling be- 



(22) 



£ ef f [p- m {tj\ = - 



Tr 

J- 1 r, 



V(t), V(t-a), 



(23) 



|0)<0|®p_ m (s)]]da, 



for this regime. The integral kernel in the above 
equation contains the two-time correlation functions, 
• (0|a m (i)oj rl (s)|0) and (0|a m (s)a^(i)|0), of the defect. 
These functions can be calculated using the dissipator 
given by Eq. ^ which, in this case, yields the inter- 
nal dynamics of the defect. By doing so and using the 
quantum regression theorem [TTJ [T^] , one finds that the 
two-time correlations decay exponentially, 

(0\a m (t)al(t - s)\0) = (0\a m (t - s)al(t)\0) = e~i s . 

(24) 

The characteristic time-scale in which the correlations 
decay is r = 2/7. In the case where the damping rate 
is large, the correlations decay rapidly, as illustrated in 
Fig. ([2]). Since the correlations are short-lived, the inter- 
actions occurs sharply at s = t. Therefore, it is safe to 
replace s in p_ m (s) in Eq. (23), by t. And by extending 
the upper limit of the integration to infinity, which con- 
cludes the Markov approximation, we can perform the 
integration with respect to s, obtaining, 



2J 2 

£eff \fi-m(t)] = — (2a m /5-m(OaL 

-aL«mP-m(i) - /5-m(t)aL«m) 



(25) 



where a m = a m _! + a TO+1 . Hence, the dynamics of the 
lossy site is no longer part of this effective description 
and the losses are described as a property of certain type 
of states, namely symmetric superpositions between the 
neighbors of the defect. 



The dissipation-free subspace in the effective picture 

Since the damping affects only symmetric superposi- 
tions, states that do not overlap with them are immune 
to particle losses. Therefore, the anti-symmetric states 
of the from, 



it 



\1>h 



|0>, 



(26) 



for any integer number n have constant population. It 
can be analytically shown that those anti-symmetric 
states identically equal to the odd modes that constitute 
the dissipation-free subspace. 
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Relation to the Quantum Zeno effect 



with 



N 



Another implication of Eq. ( 25 ) comes from the struc- 
ture of the effective dumping rate 2 J 2 j^. As the damping 
rate 7 increases, the actual losses decrease. This behav- 
ior is a manifestation of the quantum Zeno effect O [T!5] . 
This result has been experimentally discovered for ultra- 
cold molecules by Syassen et al. [2U] and is predicted 
in[Hl [TO], see also [2TJ Although this effect is not 

directly obvious from the form of the dissipator given 
by Eq. (fO), the related phenomenon is of course not a 



consequence of the adiabatic elimination. 



and 
S(i,j) 



niri \ . I nirj 
sm I sin 



x ' n— 1 v 

/ mtl \ ( mrf 

x sm sm — 

\N + 1J \N + 1 



N + l 



N- 



N 



1 ^ 

n=l 



N + l 



WKJ 

N + l 



(28) 



(29) 



IV. INTERACTIONS AND ASYMMETRY: TWO 
FACTORS THAT INDUCE LOSSES TO THE 
DISSIPATION-FREE SUBSPACE 



All the above discussion was limited to the case of a 
symmetric structure (defect at the center of the lattice) 
and no interactions. In this section we investigate how 
deviations from this ideal case can affect the protected 
modes. First we shall consider the case of a symmetric 
structure but with small yet finite interactions. Then we 
shall check the dynamics when the defect is not located at 
the center of the lattice but rather at a random position 
along it. 



Small interactions 



Phenomena that depend on interference are typically 
affected by the presence of interactions, since the latter 
destroy the phase coherence. This happens because the 
term in the Bose-Hubbard Hamiltonian ([3]), which is re- 
sponsible for the interactions, is non-linear and gives rise 
to a non-linear term in the corresponding Schrodinger 
equation. However, in the case of small enough interac- 
tions, one expects the system to behave similarly to the 
non-interacting case, at least for sort times. In order to 
examine to what extent interactions affect our findings 
of the previous sections, we compare the strength of the 
interactions to the damping rate that appears in Eq. ( 13 ) 
for an ID lattice. 



For a one-dimensional Bose-Hubbard system, the part 
of the Hamiltonian ([3| that describes the on-site particle 
interaction, in terms of the mode operators defined in 
Eq. ([©J) , reads, 



Hi, 



j,l,f)b\b p b\b } 



U 



-^J2s(i,j)bibj, 

(27) 



One can show that function S is identically equal to a 
delta function, S{i,j) = 5ij. This means that the second 



term of the right-hand side of Eq. (27) is diagonal and 



generates no hopping of particles between protected and 
unprotected modes. On the contrary, the function T is 
not a delta function and therefore the non-linear term of 



the Hamiltonian (27 1 is responsible for hopping between 



the two types of modes. 

For small enough interactions, the scattering of par- 
ticles from the protected into the unprotected modes 
occurs a lot slower than the subsequent loss of parti- 
cles from the unprotected modes. Hence, one expects 
a that particles in the protected modes still survive a lot 
longer than the other particles. To find a condition for 
when this scenario occurs, we consider the ratio between 
the scattering rate UT/2 and the effective damping rate 
r = 2j/(N + 1) as given in Eq. ( 13 ). The scattering from 



protected to unprotected modes needs to be slower than 
the losses for all modes involved. We therefore calculate 
when, 



U 
2T 



Max[|T|] < 1. 



(30) 



To this end, we numerically calculate the maximum value 

of T'(i,j,l,f) = T(i,j,l,f) as a function of the 

size of the system, N, see Fig. Q, and observe that 
max[T'] 1=3 N/2. Therefore, the ratio of the maximum 
value of the scattering over the effective damping is, 



U 



U „ r ,__ n UN 

2T Max[|r| ^2^vTi)^7 



(31) 



We thus conclude that the slow-down takes place in the 
regime where 7 3> U . An analogous treatment shows 
that this condition is also valid in two dimensions. 



B. Deviations from the symmetry 

Up to now, we have assumed that particle losses occur 
at the center of the chain, see Eq. Now, we shall place 
the losses not at the central site m but rather to a random 
position along the lattice. Instead of subscript m, we use 
r which stands for " random" . From the structure of the 
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site number 2r — 1, we define the modes, 
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FIG. 4: (Color online) The maximum value of T' as a func- 
tion of the size of the system - the total number of lattice 
sites N. The plot reveals the linear character of the function. 
Specifically, Max[T'] ~ N/2. 




FIG. 5: (Color online) Illustration of an asymmetric setup. 
The defect is located at a random position along the lattice 
rather than being at the center. For the description of the 
dynamics, we split the lattice into two parts, - each part is 
described by a different set of modes - the symmetric one that 
consists of the sites around the defect, obtaining in this way 
the structure of the symmetric setup, and the one consisting 
of the left-over sites. 



normal modes in equation Q one can see that for each 
location r, those modes with vanishing amplitude at r 
are protected against losses affecting the site r only. 

A general tendency to what degree the dissipation free 
subspace is degraded by deviations from the symmetry 
of the setup can be obtained from the following consid- 
eration. We assume, without loss of generality, that the 
defect is located at a site on the left of the center of the 
lattice. We can split the lattice into two parts, the sym- 
metric part around the defect and the remaining part 
that consists of the sites that make the right side of the 
lattice larger - an illustration is given in Fig. 

In this way, we can again retain modes that are 
symmetry-related. Specifically, for the symmetric part 
of the lattice, namely for the first site from the left to 



1 

1 



2r-l 

£ 

i=i 

2r-l 



sin 



v n— 1 



irnl 

( iml 
[~2r 



(32) 



Accordingly, for the left-over sites, namely for site 2r to 
site N, wc define the additional modes, 



VN - 2r 



N 

N 
n—2r 



il 



N -2r 

nnl 
N -2r 



ci 



(33) 



Since the lossy site r belongs to the symmetric part of the 
lattice the ladder operators a\ and a r can be expanded 
in the mode operators defined in Eq. (32 1. Therefore, the 
dissipator given by Eq. (|2j), reads in this case, 

2r-l 



C [ps(t)\ £ Sln ( T ) Sin 



{2bip s (t)b 



p s (t)b)b 



(34) 

Similarly to the ideal case, we easily observe that the 
above dissipator affects only the modes which correspond 
to odd quasi-momentum labeling number. That means 
that the modes with even quasi-momentum labeling num- 
ber that correspond to the symmetric part and all the 
modes that correspond to the left-over part are not di- 
rectly affected by the particle losses. 

Since we consider again a non-interacting case, the uni- 
tary part of the evolution of the system is governed by 
the Hamiltonian given by Eq. ([5|. In order to write this 



Hamiltonian in terms of the mode operators ( 32 ) and 
(33), we split it into three parts, 



N-l 



h.C.) - J ^ (ffl^On+l + h.c.) 



j (a: 



at 

2r- 



\0-2r 



h.C. 



(35) 

The first part of the above Hamiltonian describes the 
symmetric part of the system, the second describes the 
left-over sites and the last one gives the tunneling be- 
tween the last site of the symmetric part and the first 
one of the left-over part. Now, in terms of the modes 
(32) and (33), the Hamiltonian reads, 



2r-l / 7\ „ „ N 

fl = — 23 cos I — ) b)bi — 2 J cos I — . . ,-. , 

\ 2r J 1 ^ \N -2r ) 3 J 

j=2r 



1 = 1 
2r-l N 



etc. 



1=1 j=2r 



(36) 
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where we have defined, 



^2r{N - 2f) 



sm 



7r(2r- 1)/ 
2~r 



sm 



n2rj 



N 



2r 
(37) 

The above Hamiltonian includes a tunneling term be- 
tween the modes of the symmetric part and the left- 
over ones. This tunneling process can transform pro- 
tected to unprotected modes and destroys the effect of 
the dissipation-free subspace. However, if the tunneling 
between both parts of the lattice occurs slowly compared 
to the rate of particle losses at the defect, then one ex- 
pects a slow-down of particle losses for the modes with 
even quasi-momentum labeling number. In such a case, 
the system approximately resembles the behavior of the 
ideal case described in the first section, at least for short 
times. We therefore analyze this case in a similar way as 
the case of small interactions. The tunneling of particle 
from modes with even quasi-momentum labeling number 
to lossy modes is slow compared to the losses T affecting 
the latter provided, 



Max[\K\] < 1. 



(38) 



Inserting r = 7/ (2r) and the values of 1Z, this ratio reads, 



(39) 



This ratio depends on the relative magnitudes of the pa- 
rameters J and 7 but also on the size of the asymmetry. 
In the case of small asymmetry, where N ~ 2r, the ratio 
is approximately written, 



and since the ratio 2r/N is smaller than one, we have a 
slow-down if 7 > J. 



CONCLUSION 



In this article we have considered the Bose-Hubbard 
model with particle losses at one lattice site. For the non- 
interacting case, we found that particles in normal modes 
with vanishing amplitude at the dissipative defect are not 
affected by the localized losses. For a one-dimensional 
model with the lossy site exactly in the center of the 
chain, half the modes thus form a dissipation free sub- 
space. This behavior can be attributed to a destructive 
interference of particles tunneling into the defect. Fur- 
thermore a fraction of the particles can propagate across 
the dissipative defect even if the rate of tunneling be- 
tween adjacent lattice sites is much slower than the loss 
rate at the defect. 

To estimate the robustness of the features we predict, 
we have analyzed the effect of small particle interactions 
and deviations from a perfectly symmetric setup. 

Our findings could be studied experimentally with ul- 
tracold bosonic atoms in an optical lattice, where an elec- 
tron beam on a single lattice site ionizes atoms that are 
then extracted by an electrostatic field. 



■ Max[\TZ\] 



7 



N. 



(40) 



This indicates that there is no slow down of particle losses 
as long as 7 ^> J\fN. In the opposite case, where the 
asymmetry is large, i.e. N > 2r, we approximately have, 



J 



Max.[\n\] 




(41) 
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